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ABSTRACT 


The  displacement  equations  of  motion  for  a  6-constant  centro- 
symmetric  cubic  crystalline  material  are  developed.  Also,  the 
determinantal  equations  for  velocities  of  plane  elastic  waves  are 
developed.  Previously  determined  elastic  constants  are  used  to 
calculate  the  phase  velocities  in  various  directions  of  a  cubic 
crystal.  Normal  mode  vibrations  of  cubes  and  infinite  prisms  are 
theoretically  considered.  Dilatational  modes  of  these  bodies  are 
determined  and  they  demonstrate  that  there  are  no  couple-stress 
effects  in  these  cases.  Normal  mode  shapes  and  frequencies  for  the 
equivoluminal  modes  are  obtained  but  without  consideration  of  couple- 
stresses. 

Free  vibration  experiments  were  conducted  on  a  two-dimensional 
model  which  was  considered  as  a  slice  of  an  infinite  prism  and  on  a 
three-dimensional  cubic  model.  Experimentally  determined  frequencies 
and  mode  shapes  are  studied  in  terms  of  the  theoretical  results. 


INTRODUCTION 


In  a  previous  study  [1]  ,  a  physical  model  of  a  so-called  6- 
t  constant  centro- symmetric  cubic  crystal  was  introduced  in  order  to 

illustrate  the  phenomenon  of  couple-stresses  and  to  determine  the 
numerical  values  of  the  elastic  constants.  As  concluded  from  that 
work  the  constants  were  satisfactorily  determined,  but  it  was  then 
considered  worthwhile  to  perform  normal  mode  vibration  experiments 
to  further  explore  the  theory  and  examine  its  completeness.  In 
addition,  it  is  thought  to  be  useful  to  investigate  the  nature  of 
elastic  wave  propagation  in  such  materials.  For  these  purposes  the 
present  research  has  been  conducted.  Consequently  the  equations  of 

* 

motion  for  the  models  will  be  developed  and  solved  for  the  cases  of 
*  interest. 


EQUATIONS  OF  MOTION 

The  three-dimensional  stress  equations  of  motion  in  terms  of 
force-stress  and  couple-stress  may  be  written  down  directly  from  the 
equations  of  equilibrium  which  were  previously  given  [1]: 


Numbers  in  brackets  refer  to  references  at  end  of  report. 
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On  account  of  nonsymmetry  of  the  force-stress  tensor  we  may  write 


it  as, 
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The  couple-stress  tensor  may  be  written  as  sum  of  deviatoric  and 
spherical  parts  as  follows: 
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where  q..  stands  for  deviatoric  part  and  repeated  index  means  summa¬ 
tion.  Now  the  constitutive  equations  for  a  six-constant  centro- 
symmetric  cubic  material  may  be  written  as  follows; 
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In  the  case  of  isotropy,  the  following  relations  hold; 
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2 

The  strains  and  curvatures  may  be  expressed  in  terms  of  displace¬ 
ment  functions  u,  v,  w  as  follows: 
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where  the  03*  s  are  the  rotations  whose  components  are: 
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On  substitution  of  equations  (4),  (5),  (6)  and  (7)  into  equation  (1) 
and  after  proper  grouping  the  displacement  equations  of  motion  are  as 
follows: 
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PROPAGATION  OF  ELASTIC  WAVES 


The  differential  equations  of  motion  (8)  can  now  be  used  to  study 
the  propagation  of  elastic  waves.  The  displacements  for  plane  waves  in 
an  infinite  elastic  media  may  be  assumed  to  be  of  the  usual  form  [2], 

9*Tf 

i  -r~  (ix  +  my  +  nz  -  Ct) 

(U,v,w)  =  (U,V,W)e  (9) 

where  U,  V,  W  are  the  components  of  displacement  amplitude,  X  is 
the  wave  length,  C  is  the  velocity  of  propagation,  (f,  m,  n)  are 
direction  cosines  of  the  normal  to  wave  front. 

On  substitution  of  the  displacement  functions  (9)  into  equations 
of  motion  (8),  we  obtain  the  following  three  equations  in  terms  of  the 
amplitude  components  U,  V,  W: 

(ru  -  pc2)u  +  r12v  +  r13w  =  o 

r21U  +  (r22  "  pC-2)V  +  r23W  =  0  (10) 

r31U  +  r32V  +  (r33  "  pc2)W  =  0 
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The  condition  for  nontrivial  solutions  of  equations  (10)  is  that  the 
determinant  of  coefficients  should  vanish, 


8 


0 


(12) 


11 


-  pc2) 

ri2  ri3 

r21 

(r22  pC  )  r23 

r 

— - - 

“3-2 - - —  ^C33  pC 

Solution  of  equation  (12)  for  an  arbitrary  direction  (i,  m,  n) 
of  the  normal  to  wave  front  and  for  nonzero  values  of  the  amplitude 
components  U,  V  and  W  is  very  laborious  and  unnecessary  for  present 
purposes.  In  order  to  investigate  the  influence  of  couple-stress  on 
the  wave  velocities  and  to  identify  the  associated  nature  of  deformation 
plane  wave  normals  may  be  chosen  in  the  directions  shown  in  Figure  1. 
These  directions  are  those  of  the  cubic  axis  [100],  the  face  diagonal 
[110]  and  the  body  diagonal  [111].  By  appropriate  selection  of  U,  V 
and  W  the  phase  velocities  can  be  determined  as  definite  functions  of 
the  elastic  constants.  The  phase  velocities  in  each  of  the  above 
mentioned  directions  are  separately  treated  for  the  following  cases. 

Case  (1):  Wave  Normal  Parallel  to  Cubic  Axis  [100] 

In  this  direction, 

&  =  1  ,  m  =  n  =  0 

If  we  take  all  the  displacement  components  as  nonzero,  then  there  will 
be  three  velocities  in  each  direction. 

However,  if  we  consider 
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then  equation  (12)  will  reduce  to  the  following: 


(ru  -  pc2)  (r22  -  pc2)  -  r22  =  o 
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from  which  there  results  the  following  relations, 
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where  and  are  the  phase  velocities.  Also,  another  form  of 

displacement  amplitudes  which  will  yield  a  simple  velocity  equation  is, 
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where  is  the  phase  velocity.  It  may  be  noticed  from  equations  (14) 

and  (15)  that  . 
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Case  (2):  Wave  Normal  Parallel  to  Face  Diagonal  [110] 


In  this  direction, 
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Case  (3);  Wave  Normal  Parallel  to  Body  Diagonal  [111] 
In  this  direction, 


£  =  m  =  n  = 
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and  for 


U  =  V  t  0 
W  t  0 


the  wave  velocities  are  given  by: 
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and  also  for 
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the  wave  velocity  is  given  by 
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Note  that  in  this  case. 


Density  of  the  model  is  assumed  to  be  the  ratio  of  total  mass  of 

the  unit-cell  to  the  volume  of  unit-cell.  The  value  for  the  model  under 

consideration  is,  p  =  - — — —  l^s.:  >  where  h  (3  inches)  is 

h  in. 

the  center  to  center  distance  between  the  blocks.  For  the  case  of 
X  =  h  ,  the  wave  velocities  are  tabulated  in  Table  I.  If  the  couple- 
stress  effect  is  neglected  all  the  velocity  expressions  reduce  to  those 
of  the  classical  case  for  a  cubic  crystal  [2], 
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NORMAL  MODE  VIBRATIONS 


Now  consideration  will  be  given  to  the  normal  mode  vibrations  of 
the  crystals.  The  vibrating  body  is  assumed  to  have  stress-free  boundary 
conditions.  Theoretical  studies  of  normal  mode  vibrations  of  an  infinite 
prism  on  the  basis  of  plane  strain  equations  and  of  a  cube  on  the  basis 
of  three-dimensional  equations  are  considered. 

a 

1.  Plane  Strain  Vibrations  of  Infinite  Prisms 

The  displacement  equations  of  motion  (8)  for  plane  strain  reduce 

to, 
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Taking  the  origin  at  the  center  of  the  model,  we  may  write  down 
the  stress-free  surface  conditions  which  should  be  satisfied  at 
x  =  +  a  and  y  =  +  a  as  follows; 
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(a)  A  Dilatational  Mode 

Since  the  dilatational  mode  consists  only  of  extension  and  con¬ 
traction  of  principal  line  elements  of  the  material,  the  frequencies 
are  not  influenced  by  couple-stress.  Displacement  functions  which  will 
satisfy  equations  (20)  and  (21)  for  the  case  of  vanishing  Poisson's 
ratio  are  found  to  be, 
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with  m  representing  the  odd  integers.  The  above  functions  will  describe 
symmetric  modal  shapes  with  respect  to  the  coordinate  axes.  Asymmetric 
modes  are  given  by 
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with  m  representing  the  even  integers. 


The  frequency  equation  corresponding  to  functions  (22)  and  (23) 
are  found  to  be. 
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where  m  =  1,  2,  3, 


(b)  An  Equivoluminal  Mode  Without  Couple-stresses 

Without  consideration  of  the  couple-stresses,  the  differential 
equations  of  motion  (20)  and  boundary  conditions  (21)  reduce  to 
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Displacement  functions  which  will 
for  symmetric  modes  are  seen  to  be  the 


satisfy  equations 
following: 


(25)  and  (26) 
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with  m  representing  the  odd  integers.  The  asymmetric  modes  are  given 
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with  m  representing  the  even  integers. 

The  frequency  equation  corresponding  to  functions  (27)  and  (28) 
are  found  to  be. 
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2.  Vibration  of  a  Cube 

All  three  displacement  functions  u,  v,  w  are  involved  in  this 
case  and  the  equations  of  motion  are  equations  (8) .  Considering  the 
center  of  the  cube  as  origin  of  the  coordinate  system,  the  conditions 
of  stress- free  surfaces  at  x  =  +  a  ,  y  =  +  a  and  z  =  +  a  are  the 
following: 

s  s  <?  1  Sqd  dqd 

2  OX  dy  bz  bz 

,  dq^  dq^  dq^  d 

.  ±  (_3  +  +  Jzl  .  32_)n  .  0 

2  by  <Dz  dy 


s  s  ,  Sq^  dqd  <5qd  ^-d 

■La  +  Tim  +  Tf  n  -  i  (-22  +  _S!2  +  -is  .  aCL, 


2  bx  dy  Bz  Bz 

,  Bqd  Bqd  v-d 

+  i  (-s  +  Jh  +  _!h  .  S_)n  . 

2  Bx  by  bz  bx 


T  i  +  T 
xz  yz 


s  c  1  dq^  Bq^  Bqd  -n— d 

T>  +  T2S2n  +  i  ( JS  +  JH  +  J2Z  .  SS_} 

2  dx  cly  bz  cty 

■  1  oq^  dq^  dq^  --d 

-  -L  (  Jh  +  .12  +  Js£  . 

2  bx  by  bz  Bx 


m  =  0 


d»d  d  — d . 

q„&  +  q  m+q  n  ■  Q  I 
xx  nyx  Xzx  x 

d-d  d  — d 

qt  +  q  m+q  n-Qm 
^xy  nyy  Mzy  x 

d  -  d  d  — d 

q„_i>  +  q  m+q  n-Qn 
xz  ^yz  ^zz  x 


where  l>,  m,  n  are  direction  cosines  of  unit  normal  to  boundary 

^  =  02  +  <yJ  +  V+  <4+ V)fc  +  <«£.+  O  +  <£  "  d) 


yz  ^zy 
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(a)  A  Dilatational  Mode  of  a  Cube 

In  the  case  of  zero  Poisson's  ratio,  displacement  functions  which 
will  satisfy  equations  (8)  and  (30)  are  the  following: 


.  .  >  mrt  x  icot 

u  =  A  sm(— — )x  e 

2a 


.  .  ,  mrt  .  icot 

v  =  A  sin( - )y  e 

2a 


.  .  .  nut  N  icot 

w  =  A  sm( - )z  e 

2a 


(31) 


with  m  as  odd  integers.  The  above  functions  will  describe  symmetric 
mode  shapes  with  respect  to  the  coordinate  planes.  Asymmetric  modes 
are  given  by, 


*  .nm_ .  icot 

u  =  A  cos( - )x  e 

2a 


.  .nut  .  icot 
v  =  A  cos( - )y  e 

2a 


.  .  mit  .  icot 

w  =  A  cos(- — )z  e 

2a 


(32) 


with  m  as  even  integers. 


The  frequency  equation  corresponding  to  functions  (31)  and  (32)  is 
as  follows: 


pco 


2 


su(— ) 

2a 


where  m  =  1,  2,  3,  ...  . 


(33) 
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As  in  the  two-dimensional  case  the  frequencies  of  the  dilat.ational 
modes  of  the  cube  are  not  influenced  by  couple-stress. 

(b)  An  Equivoluminal  Mode  of  a  Cube 

Without  considering  the  effect  of  couple-stress,  the  equations  of 
motion  (8)  reduce  to  the  following: 


-v2  >2  v2  S..  >2  .2  >2  n2  n2 

o  u  .  ,ov  .  o  w  ^  44  ,o  u  o  u  o  v  .  0  w  ,  du 

S  — n  +  S.  9  ( -  +  - )  +  -  (  r  +  r  + - + - )  =  p  — r 

dx  '  dxdy  <5xbz  2  dy  dz  dxdy  ()xc)z  dt 


■\2  n2  \2  S.,  n2  n2  n2  s2  n2 

ov.„  ,  a  w  .  a  u  v  .  44  ,o  v  .  o  w  .  o  v  ,  o  u  x  ov 

s,  ■,  ~7+s19  ( - +  — )  +  —  (-7  + - +  —o  + - )  =  p  — 7 

dy  dyBz  dydx  2  ciz  c)zdy  c)x  dxdy  Bt 


>.2  n2  S.  .  ^2  -\2  v.2  -v2 

ou  ,  ovN,  44/Ow,  ou  ,  o  w  ,  ov 


S1  1  - 9  +  S1  9  ( - +  - )  +  -  (  9  +  -  +  9  +  - ) 

dz  dzdx  dzdy  2  Bx  dzck  dy  dzcfy 


and  boundary  conditions  (30)  reduce  to: 


the  following: 


u  = 

sin(-^-)x[A„ 

2a 

,mx  \  * 

cos ( - )y  -  A 

2a  1 

,  nut .  ,  icut 
cos( - )z]e 

2  a 

V  = 

sin(— )y[Ax 
2a  L 

cos(-^-)z  -  A0 
2a  J 

.mrt  x  ■,  icot 
cos( - )x]e 

2a 

(36) 

w  = 

sin(— )z[A„ 

cos(-^-)x  -  A. 

.  mx  .  .  i(ot 
cos( - )y]  e 

2a 

2a 

2a 

with  m  as  odd  integers.  Asymmetric  modes  are  given  by 


cos(-^-)x[A~ 

sin(— )y  -  A2 

.  .  mx  s  .  icot 
sm  ( - )z]  e 

2a  J 

2a 

2a 

cos(— )y[A1 

sin(-^~)z  -  A- 

.  ,mx  N  ,  i(JDt 
sm( - )x]  e 

2a 

2a 

2a 

cos(-^^-)z[  A_ 

sin(-^-)x  -  A. 

.  ,  mx  N  n  icot 
sm( - )y]e 

2a 

2a 

2a 

with  m  as  even  integers. 


(37) 


Frequency  equation  corresponding  to  equations  (36)  and  (37)  is 
seen  to  be : 

9  r  2 

P»  =  <SU  -  S12)(— )  (38) 

2  a 

where  m  =  1,  1,  3}  ...  . 
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EXPERIMENTAL  STUDY  OF  FREE  VIBRATIONS 


Using  the  experimentally  determined  elastic  constants,  one  can  use 
the  solutions  of  the  equations  of  motion  to  predict  the  normal  mode 
shapes  and  corresponding  frequencies  of  the  models.  It  is  thought  to 
be  of  considerable  value  for  the  development  of  the  theory  of  couple- 
stresses  to  have  knowledge  of  the  actual  vibration  characteristics  of 
the  physical  models.  The  experiments  and  results  will  now  be  described 
in  some  detail. 

1.  Experimental  Method 

The  experimental  models  used  in  the  vibration  experiments  were  sus¬ 
pended  by  flexible  strings  from  rigid  steel  frames  which  were  bolted  to 
a  massive  reinforced  concrete  base  3  X  3  X  10  feet.  The  freely  suspended 
models  will  give  the  stress-free  surface  conditions  for  free  vibration 
of  solids.  Mounting  of  models  and  necessary  electrical  connections  to 
excite  the  models  are  shown  in  Figures  2  and  6. 

Experimental  method  that  has  been  used  is  similar  to  the  one  dis¬ 
cussed  in  a  paper  by  Hoppmann  [3],  Normal  mode  vibrations  were  excited 
in  the  solids  by  an  electromagnet  mounted  solidly  at  the  appropriate 
location.  Since  the  experimental  models  are  partly  made  out  of  aluminum 
blocks,  it  was  necessary  to  attach  a  laminated  iron  piece  to  the  model 
at  a  suitable  point  to  facilitate  the  magnetic  excitation.  The  response 
of  the  model  was  picked  up  by  crystal-type  phonograph  pick-up. 

The  magnet  is  excited  from  an  audio-oscillator  which  feeds 
through  a  resonant  type  RC  circuit  in  order  to  increase  the  driving 
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force.  On  account  of  the  nature  of  the  magnet,  the  frequency  from  the 
oscillator  is  one-half  the  frequency  of  the  driving  force  of  the  magnet. 
The  output  from  the  oscillator  is  connected  to  one  set  of  plates  of  a 
cathode  ray  oscilloscope  and  the  magnetic  drive  is  connected  across  the 
other  set  of  plates.  Since  the  frequency  ratio  is  2  to  1,  a  figure 
eight  Lissajous  is  shown  on  the  screen  of  the  scope.  Crossing  a  nodal 
surface  with  the  pick-up  turns  the  figure  eight  Lissajous  into  a  horse¬ 
shoe  shaped  figure  on  one  side  of  the  surface  and  inverts  the  figure  on 
the  opposite  side.  So  by  moving  the  crystal- type  pick-up  over  the 
models  the  nodal  surfaces  and  hence  the  nodal  patterns  were  determined. 
This  is  a  precise  method  for  the  determination  of  the  modal  shapes.  It 
is  important  to  determine  the  mode  shapes  and  their  corresponding  fre¬ 
quencies  simultaneously,  especially  since  the  spectrum  of  frequencies 
may  be  closely  packed  for  any  elastic  body  of  the  nature  under  investiga¬ 
tion. 

2.  Experiments  on  Model  in  Shape  of  a  Slice  of  an  Infinite  Prism 

The  two-dimensional  model  was  considered  as  a  slice  of  an  infinite 
prism.  The  model  was  freely  suspended  (Figure  2)  from  the  solid  frame 
by  means  of  two  flexible  s -rings  so  that  no  constraint  was  introduced. 

The  model  was  excited  so  that  motion  out  of  the  plane  of  the  model  was 
prevented.  The  vibrations  are  then  approximately  those  of  plane  strain. 

A  schematic  diagram  of  the  experimental  apparatus  with  some  details 
of  model  and  electrical  connections  is  shown  in  Figure  2.  The  magnet 
was  mounted  on  the  solid  frame  and  was  easily  movable  in  any  direction. 
The  model  was  excited  at  different  points  and  the  mode  shapes  and  fre¬ 
quencies  were  determined. 
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After  determining  the  mode  shapes,  an  attempt  was  made  to  classify 
them.  It  is  well-known  tnat  the  complete  classification  of  normal  mode 
vibrations  for  elastic  solids  has  never  been  accomplished.  Some  attempts 
have  been  made  with  some  degree  of  success  as  discussed  by  Ekstein  and 
Schiffman  [4], 

On  comparing  the  experimentally  determined  mode  shapes  with  the 
theoretical  shape  functions  it  is  considered  that  they  represent  the 
dilatational,  equivoluminal  and  face-shear  types.  In  the  classification 
of  Ekstein  and  Schiffman  [4],  the  fundamental  modes  of  dilatational, 
equivoluminal  and  face- shear  type  are  referred  to  as  the  breathing, 
longitudinal  and  shear  modes,  respectively. 

The  face-shear  type  modes  were  the  lowest  in  frequency.  The  fre¬ 
quencies  of  the  equivoluminal  and  dilatational  modes  were  rather  close 
to  each  other.  The  mode  shapes  along  with  the  corresponding  fi-equencies 
are  shown  in  Figures  3,  4  and  5.  It  was  found  that  the  experimental 
nodal  lines  shown  in  Figure  3  were  approximately  the  same  as  that  given 
by  shape  functions, 

u  =  A  sin(-^—)y 
2a 

v  =  A  sin(-^-)x 
2a 

for  odd  values  of  m  and 
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for  even  values  of  m  .  Accordingly,  the  mode  shapes  are  sketched  from 

these  functions.  Calculated  and  experimental  frequencies  are  tabulated 

.  ,  .401  X  10  4  lbs. sec. ^ 

m  Table  II.  The  average  density,  p  =  -  - r - ,  of  the 

a  h  in.^ 

model  was  used  in  calculation. 

3.  Experiments  on  Three-Dimensional  Model 

The  model  was  freely  suspended  from  a  solid  frame  by  means  of  three 
inextensible  strings.  The  method  of  exciting  the  model  was  similar  to 
that  discussed  previously.  A  photograph  of  the  apparatus  is  shown  in 
Figure  6. 

The  frequencies  and  mode  shapes  were  determined  by  exciting  the 
model  at  different  points.  An  attempt  has  been  made  to  classify  the 
mode  shapes.  Besides  the  analogous  shapes  of  two-dimensional  models 
there  were  certain  modes  which  were  not  identified.  The  classified 
mode  shapes  and  frequencies  are  shown  in  Figures  7,  8  and  9.  The  un¬ 
classified  mode  shapes  and  frequencies  are  given  in  Figure  10.  Theo¬ 
retical  frequencies  are  calculated  on  the  basis  of  average  density  of 

~4  2 

the  model,  p  =  - s  .  Calculated  and  experimental 

a  h  in.4 

frequencies  are  shown  in  Table  III. 


DISCUSSION  AND  CONCLUSION 


For  the  models  which  were  studied,  it  has  been  found  that  the  wave 
velocity  shows  dispersive  effects  because  of  couple-stresses  and  as  the 
wave  length  becomes  indefinitely  decreased  the  phase  velocity  becomes 
unbounded.  Also,  it  was  found  that  the  velocity  of  a  purely  dilatational 
wave  is  uninfluenced  by  the  couple  rigidity  characteristic  of  the  model. 

As  is  well-known,  the  classification  of  normal  mode  vibrations  for 
elastic  solids  is  incomplete.  Ekstein  and  Schiffman  [4]  have  attempted 
to  classify  the  modes  of  isotropic  cubes  by  using  group  theory,  but  only 
with  a  modest  degree  of  success.  Classification  of  the  modes  for  centro- 
symmetric  cubic  crystals  with  or  without  couple-stress  has  never  been 
attempted  so  far  as  is  known.  In  the  present  paper  some  attempt  was 
made  to  classify  the  modes  on  the  basis  of  experimental  and  theoretical 
results. 

As  would  be  anticipated  from  analysis  of  normal  mode  /ibration  the 
dilatational  modes  are  unaffected  by  the  couple  rigidity  factor.  Equi- 
voluminal  and  face-shear  type  modes  do  show  couple-stress  effects. 

Since  the  experimental  model  is  inhomogeneous  on  account  of  the 
manner  in  which  its  structured  form  and  mass  distribution  have  been  de¬ 
signed,  the  vibration  experiments,  in  relation  to  the  theory,  can  only 
be  expected  to  be  meaningful  for  wave  lengths  which  are  longer  than  the 
size  of  a  typical  cell  unit.  Also,  for  the  higher  modes  it  should  be 
expected  that  there  will  be  some  influence  of  rotatory  inertia  of  the 
discrete  blocks  which  does  not  arise  in  the  theory  of  continuous  media. 
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From  the  vibration  experiments  which  were  conducted,  it  was  found 
that  the  face-shear  type  modes  fall  into  a  lower  set  of  frequencies, 
the  highest  of  which  is  about  250  cps,  and  the  dilatational  and  equi- 
voluminal  modes  are  in  a  higher  frequency  set  which  starts  at  about 
1100  cps.  Also,  it  was  seen  that  the  experimentally  determined  fre¬ 
quencies  of  the  lower  modes  check  the.  theoretical  values  to  within  5%. 

As  expected,  the  frequencies  do  not  check  that  well  for  the  higher  modes 
for  the  reason  given  above. 

Unfortunately,  since  it  was  not  possible  to  solve  the  equations  of 
motion  for  those  free  vibrations  of  the  models  which  would  be  expected 
to  show  couple-stress  effects,  it  was  then  not  possible  to  predict  the 
influence  of  couple-stress  on  the  frequencies  of  vibration.  However,  it 
is  interesting  to  observe  that  there  was  close  agreement  between  experi¬ 
mentally  determined  modes  of  vibration  and  the  theoretical  modes  of 
vibration  which  obviously  do  not  involve  couple-stresses. 
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TABLE  I 


PHASE  VELOCITIES  OF  WAVE  PROPAGATION  IN  VARIOUS 
DIRECTIONS  OF  THE  PHYSICAL  MODEL 


Displacement 


[100]  Direction:  i 


UjEO,  V  /  0,  W  =  0 


11 


=  n  =  0 


U  =  0,  V  =  0,  WM 


Velocity  _  _1_  l_  j44  ,  Ci it2. 

C2  12  p  C  2  +  Vr,) 

S44  ^  °-’  C  =  °  S44  *°>  C  ^  ° 

X  =  3" 

471  ft. /sec.  4270  ft. /sec. 


[110]  Direction:  l  =  , 


Displacement  U  ^  0,  V  /  0,  W  =  0 


c  I  I  (Sn  *  S44a 

C1  "  12  p  ^  2  ; 

=  3780  ft. /sec. 


c  .  i  i.  (!« + 

c3  -  n  ,  p  (  2  +  ir> 


C3  =  C2 


n  =  0 


U  =  0,  V  =  0,  W  /  0 


C3  =  " 


c  _  J_  _L  (S  +  iCn  ) 

C9  "  io  >  ^11  +  —J-) 


c  ?  2 

1  °44  ,  Crt  ,  3A  it. 

—  (. - +  x  +  —X) 

p  2  2XZ  4  Xz 


2  12  1  2p  H  X 


C  =  0 


C  /  0 

su  *  0 

C  =  A  =  0 
S44  *  ° 

C  jt  0,  A  =  0 
S44*° 

C  /  0,  A  /  0 

S44  *  ° 

X  =  3" 

5680  ft. /sec. 

471  ft. /sec. 

X  =  3" 

3050  ft. /sec. 

X  =  3" 

3060  ft. /sec. 

Displacement 


[111]  Direction:  Z  =  m  =  n  = 


U  =  V  /  0>  W  /  0 


_  1  f  1  .fll  ,  _2_  c  ,  J_  A*\ 

1  12  p  3  3  44  3  XZ 


S^i^O,  S^/0,  A-0  S^/0,  S^?£0jA/0 


X  =  3" 

3120  ft. /sec. 


U  =  -V,  W  =  0 


.  ,S11  ,  S44  ,  2C3t2  ,  5 Ait 

”*  (  T  I  11  ft  »  r 

>  3  6  9)/  6X‘ 


*11*’**#’  sn*°’ s44|i0,  sii,l0>Vi0> 

C=A=0  C/0,  A=0  C/0,  A=0 


X  =  3" 

3080  ft. /sec  3680ft. /sec. 


TABLE  II 


FREQUENCIES  OF  TWO-DIMENSIONAL  MODEL 


(a)  A  Dilatationel  Mode 


m 


Frequency  (cps) 


1 


2 


3 


4 


Experimental 


1150 


1880 


2900 


3600 


Theoretical 


1195 


2390 


3585 


4780 


(b)  An  Equivoluminal  Mode 


m 

Frequency  (cps) 

Experimental 

Theoretical 

TABLE  III 


FREQUENCIES  OF  THREE-DIMENSIONAL  MODEL 


(a)  A  Dilatational  Mode 


m 

Frequency  (cps) 

Experimental 

Theoretical 

n 

1520 

1610 

2 

2800 

3220 

(b)  An  Equivoluminal  Mode 


m 

Frequency  (cps) 

Experimental 

Theoretical 

1460 

1610 

2 

2640 

3220 

FIGURE  1.  Various  Directions  of  Wave  Normal 


Considered  in  the  Calculation  of  Phase  Velocities 


-magnetic 


Figure  2.  APPARATUS  FOR  VIBRATION  OF  TWO-DIMENSIONAL  MODEL 


MODE  SHAPES  AND  FREQUENCI 
TWO  DIMENSIONAL  MODEL 
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FREQUENCE  OF  TWO  DIMENSIONAL  MODEI. 


Figure  6.  APPARATUS  FOR  VIBRATION  OF  THREE-DIMENSIONAL  MODEL 
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Figure  LO.  NODAL  POINTS  AND  FREQUENCIES  OF  THREE-DIMENSIONAL  MODEL 
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